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Recent experimental studies of Bose-Einstein Correlations in Z frag- 
mentation are reviewed in view of the need to understand their appar- 
ent suppression for pions originating from different W's. Particular fea- 
tures discussed are source elongation, position-momentum correlation, non- 
Gaussian shape of the correlator, transverse-mass dependence, density de- 
pendence and dilution, space-time shape of the emission function, neutral- 
pion and genuine higher-order correlations. 



1. Introduction 

As proposed by Hanbury Brown and Twiss [1] in 1954, the (angular) 
diameter of stars and radio sources in the universe was successfully deter- 
mined by measuring the intensity correlations between separated telescopes. 
Likewise, in particle physics, one can in principle use Bose-Einstein corre- 
lations between identical bosons to measure the space-time structure of the 
region from which particles originate in a high-energy collision [2] , provided 
these bosons are produced incoherently. 

The first experimental evidence for Bose-Einstein correlations in particle 
physics dates back to 1959 when, in pp annihilation at 1.05 GeV/c, Gold- 
haber et al. [3] observed an enhancement at small relative angles in like-sign 
pion pairs not present for unlike-sign pairs. More recently, Bose-Einstein 
correlations have been exploited in hadron-hadron, hadron-nucleus, nucleus- 
nucleus, e + e~ and lepton-hadron collisions to obtain surprisingly detailed 
information on the space-time development of particle production. 

The recent revival of interest comes from various directions: 
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1. Their application to determine the space-time development of a particle 
collision. 

2. Their influence on the measurement of effective masses, in particular of 
the W mass at LEP2 [4, 5]. 

3. Their role in the phenomenon of intermittency [6, 7]. 

4. Their possible effect on multiplicity distribution and single-particle spec- 
tra [8, 9, 10]. 

In this paper, we shall review recent experimental studies on the first 
point, in particular for e + e~ collisions leading to hadronic final states at 
the Z energy. We consider this information crucial for an understanding of 
the underlying dynamics, and in particular of the apparent suppression of 
Bose-Einstein correlations of pions originating from different W's within the 
same event. 

2. The correlation formalism 

We start by defining symmetrized inclusive ^-particle distributions 

1 da g (p 1 ,...,p g ) 
Pq{Pi, ■ ■ ■ ,Pq) = - q , (1) 

1 

where o~ q (pi, . . . ,p q ) is the inclusive cross section for q particles to be at 
pi,...,p q , irrespective of the presence and location of any further particles, 
Pi is the (four-) momentum of particle i and <7tot is the total hadronic cross 
section of the collision under study. 

For the case of identical particles, integration over an interval Q in p- 
space yields 

Pi{p)dp = (n) , 

P2(pi,P2)dpidp 2 = {n(n - 1)) , 

/ dpi... I dpgpgipi, ...,p q ) = (n(n - 1) . . . (n - q + 1)) , (2) 
Jn Jn 

where n is the multiplicity of identical particles within 17 in a given event 
and the angular brackets imply the average over the event ensemble. 

Besides the interparticle correlations we are looking for, the inclusive q- 
particle number densities p q (pi, ■ ■ ■ ,p q ) in general contain "trivial" contribu- 
tions from lower-order densities. It is, therefore, advantageous to consider a 
new sequence of functions C q (pi, . . . ,p q ) as those statistical quantities which 
vanish whenever one of their arguments becomes statistically independent 
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of the others. Deviations of these functions from zero shall be addressed as 
genuine correlations. 

The quantities with the desired properties are the correlation functions 
- also called (factorial) cumulant functions - or, in integrated form, Thiele's 
semi-invariants. [11] A formal proof of this property was given by Kubo.[12] 
The cumulant correlation functions are defined as in the cluster expansion 
familiar from statistical mechanics via the sequence: [13, 14, 15] 

Pi(l) = C'i(l), (3) 
p 2 (l,2) = Ci(l)Ci(2) + C 2 (l,2), (4) 
p 3 (l,2,3) = Ci(l)Ci(2)Ci(3) + Ci(l)C 2 (2,3) + Ci(2)C 2 (l,3) + 

+ Ci(3)C 2 (l,2)+C 3 (l,2,3); (5) 

and, in general, by 

p m (l,...,m) = £ E [CiO- - -C 1 ()][C 2 U---C 2 (,)] --- 
ft} m P erm - factors h factors 

•••[C m (,...,)---C m (,. ..,)]. (6) 

V v ' 

factors 

Here, U is either zero or a positive integer and the sets of integers {k} m 
satisfy the condition 

m 

'Y^il i =m. (7) 

i=\ 

The arguments in the Q functions are to be filled by the m possible momenta 
in any order. In the above relations we have abbreviated C q (pi, . . . ,p g ) to 
C q (l, 2, . . . , q); the summations indicate that all possible permutations must 
be taken (the number under the summation sign indicates the number of 
terms). The sum over permutations is a sum over all distinct ways of filling 
these arguments. For any given factor product there are precisely [14] 

m! (8) 



[(iiyi(2\y*---(m\y™]i 1 u 2 \...i rn \ 

terms. 

The relations (6) may be inverted with the result: 

C 2 (l,2) = /> 2 (l,2)- Pl (l) Pl (2) , 
C 3 (l,2,3) = #,(1,2,3) -^ Pl (l)p 2 (2,3)+2p 1 (l)p 1 (2)p 1 (3) , 

(3) 
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C 4 (l,2,3,4) = p4(l,2,3,4)-^ Pl (l)p3(l,2,3)-^p 2 (l,2)p 2 (3,4) 

(4) (3) 

+ 2j> 1 (l)p 1 (2)p 2 (3,4) - 6^i(l)pi(2)^i(3)pi(4), (9) 

(6) 

etc. Expressions for higher orders can be derived from the related formulae 
given in [16]. 

It is often convenient to divide the functions p q and C q by the product 
of one-particle densities, which leads to the definition of the normalized 
inclusive densities and correlations: 

Rq(Pl,...,p q ) = p q (p q ,...,P q )/pi(pi) ...pi{p q ), (10) 
K q (pi,...,p q ) = C q (p 1 ,...,p q )/p 1 (p 1 )...p 1 (p q ). (11) 

In terms of these functions, correlations have been studied extensively for 
q = 2. Results also exist for q = 3, but usually the statistics (i.e. number of 
events available for analysis) are too small to isolate genuine correlations. 
To be able to do that for q > 3, one must apply moments defined via the 
integrals in Eq.(2), but in limited phase-space cells [6]. 

3. Alternative views 

3.1. Pion interferometry 

In Fig. 1 we illustrate the production of two identical pions with mo- 
menta p 1 and p 2 , arising from two sources A and B with coordinates xa, 
xb (see also [17]). The pion wave functions can be written as 

V>ia = e ~ ip i-* A+ia , V>2B = e -ip »- XB+i/J (12) 

if -n^pjj is emitted from source A and 7r(p 2 ) from source B, where a and (3 
are arbitrary phases of the sources. 
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Figure 1. Emission of two identical bosons with momenta p 1; p 2 from two sources 
A, B. 

If vr(p 1 ) is emitted from source B and 7r(p 2 ) from source A, then indices 
A and B (and the phases) should be interchanged in (12). Since the two pi- 
ons are identical bosons and the observer cannot decide from which source a 
particular pion was emitted, the coincidence amplitude for simultaneous ob- 
servation of two pions with momenta p : and p 2 has to be Bose-symmetrized: 

Abe = ^iaY^b + V'ibV^a • (13) 
The corresponding coincidence rate is 

Ibe = \A B e\ 2 = 2 + 2 cos(Ap • Ax) , (14) 

where Ap = p 1 — p 2 and Ax = xa — xb- Note that the arbitrary phases 
a, (3 have dropped out from (14), which is valid for completely incoherent 
emission. We define as two-particle Bose-Einstein ratio i? 2 the ratio be- 
tween /be and the rate Iq which would be observed if there were no BE 
interference: 

R2 = /be//o = 1 + cos(Ap • Ax) . (15) 

From (15) it follows that i? 2 reaches a maximum value of 2 for Ap = 0. 
Furthermore, it can be seen that the momentum difference Ap probes the 
source dimensions in a direction parallel to Ap. 

We shall, however, see later on that such a simple picture cannot be 
maintained due to correlation between x and p observed in the data and 
expected from hydrodynamical models as well as from string models. 

One step more realistic than the binary source considered in Fig. 1 is a 
source with a spherically symmetric Gaussian density distribution of emit- 
ting centres [3] 

p(r) oc exp [-r 2 /(2r 2 )] , (16) 
which yields as Bose-Einstein ratio 

R 2 = 1 + exp [-r 2 Ap 2 ] . (17) 

or, in its Lorentz-invariant form, 

R 2 (Q 2 ) = 1 + eM-r 2 G Q 2 ) 
with (18) 
Q 2 =-{pi- P 2? = M 2 -Ami , 

where M is the invariant mass of the pion pair. This corresponds to a 
Gaussian shape of the source in the centre-of-mass system of the pair, where 
q = AE = 0. 
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3.2. Emission function and Wigner function 

The picture presented in Sub-Sect. 3.1 corresponds to the one-dimensio- 
nal treatment of a spherically symmetric static source. However, a high- 
energy collision is neither spherically symmetric nor static. A formalism 
particularly handy for the fully-dimensional treatment of a dynamical emit- 
ter is the so-called Wigner-function formalism [18, 19]. This is based on the 
emission function S(x,p), a covariant Wigner-transform of the source den- 
sity matrix. S(x,p) can be interpreted as a quantum-mechanical analogue 
of the classical probability that a boson is produced at a given space-time 
point x = (t, r) with a given momentum-energy p = (E, p). 

In the general case, the normalized two-particle density R 2 (l, 2) or cor- 
relation function K 2 (l,2) depend on the momentum components of par- 
ticles 1 and 2. For the study of correlations, it is conveniant to decom- 
pose the two single-particle four-vectors p\ and p2 into the average K = 
[(E\ + E 2 )/2, K = p l + p 2 )/2] and the relative momentum Q = (AE = 
Ei - E 2 ,Q = Pi - p 2 ). 

Starting from the space-time x and momentum-energy K dependent 
pion-emission function S(x, K), the normalized density in momentum space 
can be written as [19] 

^"^'^y-^'^i'- (19) 

In a Gaussian approximation around the mean space-time production point 

x, 

R 2 (Q, K) = 1 + exp[-Q M Q^((x - x)^{x - z)„)(K)] + SR 2 (Q, K) . (20) 

The variances ((x — x)^(x — x) v ) give the size of the space-time region from 
which pions of similar momentum are emitted (which, for Gaussian sources, 
coincides with the more general concept [20] of lengths of homogeneity) and 
5R 2 contains all non-Gaussian contributions, usually assumed to be small. 

Since the four-momenta pi of the two particles are on-shell, Q and K are 
in general off-shell but obey the orthogonality and mass-shell constraints 

Q ■ K = , K 2 - Q 2 /A = m 2 , (21) 

so that only 6 linear combinations of the variances are measurable [21]. If the 
source is azimuthally symmetric in coordinate space, a reflection symmetry 
is present in momentum space with respect to the plane spanned by K and 
the event axis. As a consequence, all mixed variances linear in the direction 
orthogonal to this plane ( "sidewards" ) must vanish and the correlator must 
be symmetric under Qsidc - * — Qsidc, so that only four linear combinations 
remain measurable! Note, however, that every one of them depends on K. 
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Figure 2. Space-time diagram for two ways to produce two identical bosons in the 
color-string picture [23] . 



3.3. String models 

Alternatively, Bose-Einstein correlations have been introduced into string 
models [22, 23, 24]. In these models, an ordering in space-time exists for the 
hadron momenta within a string. Bosons close in phase space are nearby 
in space-time and the length scale measured by Bose-Einstein correlations 
is not the full length of the string, but the distance in boson-production 
points for which the momentum distributions still overlap. 

Fig. 2a illustrates the production of (identical) particles 1 and 2 from a 
color string in x and t. The color field breaks up into quark- antiquark pairs 
and adjacent quarks and antiquarks recombine into mesons. The production 
of the same final state, but with particles 1 and 2 exchanged is described in 
Fig. 2b. 

In a color-string model, the (non-normalized) probability dr n to produce 
an n-particle state {pj}, j = 1, . . . n of distinguishable particles is 

dr„ = [U^NdpjS^ - rr%)\8(Vpj - P) exp(-L4 n ) , (22) 

where the exponential factor can be interpreted as the square of a matrix 
element 

M n = exp(i£A n ) , Re(0 = re , Im(£) = 6/2 , (23) 

and the remaining terms describe longitudinal phase space, with P being the 
total energy-momentum of the state. N is related to the mean multiplicity 
and b to the correlation length in rapidity. A n corresponds to the total 
space-time area covered by the color field (Fig. 2), or to an equivalent area 
in energy-momentum space divided by the square of the string tension re = 1 
GeV/fm [23]. 

The production of two identical bosons (1,2) is governed by the sym- 
metric matrix element 

M = —={M 12 + M 21 ) = — [exp(i£Ai 2 ) + exp(ifA 2 i)] . (24) 
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From Fig. 2 it is clear that there is an area difference and, consequently, a 
phase difference between M\ 2 and M 2 \ of 

AA = \A 12 -A 21 \ = ^\ Pl E 2 - p 2 Ei + ( Pl - p 2 )Ej - {E 1 - E 2 ) Pl \ , (25) 

where the indices 1,2 and I represent particles 1, 2 and system I, respectively. 
Using this matrix element, one obtains 

R BE « 1 + (cos(KAA)/cosh(bAA/2)) , (26) 

where the average runs over all I. In the limit Q 2 = —{pi — p 2 ) 2 = 0, (25) 
gives AA = and (26) Rbe = 2, in agreement with the results from the 
conventional interpretation for completely incoherent sources. However, for 
Q 2 7^ follows an additional dependence on the momentum p\ of the system 
I produced between the two bosons. 

Corrections to (26) are necessary due to non-zero mass and transverse 
momentum of quarks and due to the contribution of resonances to the pro- 
duction of particles of type 1, 2. 

The model can account well for most features of the e + e~ data [25, 26, 
27], including the approximately spherical shape of the BE effect. More 
recently, the symmetrization has been generalized to more than 2 identical 
particles [28]. 

4. Recent experimental results 

4-1. Existence 

Bose-Einstein correlations are by now a well established effect in the 
hadronic final states of Z decay [29, 30, 31]. A clear enhancement is observed 
in R 2 at small Q. This is not a trivial observation since, according to the 
pion interferometry interpretation of Sect. 3.1, this would require at least 
partially chaotic pion production. 

If present in hadronic Z final states, there is no reason to expect it to be 
absent in hadronic W final states (intra-W BEC), and a signal consistent 
with that of Z into light quarks is indeed established [32] . Examples for R 2 
as a function of Q are given in Fig. 3, both for W and Z fragmentation. 

The important question is that of BEC between pions each originating 
from a different W in fully hadronic W + W~ final states (inter-W BEC). If 
existent, such a correlation would, on the one hand, cause a potential bias 
in the mass determination of the W [4, 5]. On the other, it could serve 
as a pion-interferometry laboratory for the measurement of the space-time 
development of W fragmentation into pions. The recent status of the search 
for inter-W BEC is covered in [33]. 
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Figure 3. The Bose-Einstein correlation function R 2 for (a) the fully-hadronic WW 
events, and (b) the semi-hadronic WW events. In (b) the full histrogram is for the 
light-quark Z decay sample and the dashed histrogram is for a sample containing 
all hadronic Z decays. Also shown are Gaussian fits to the WW data [32]. 

For a detailed understanding of W + W~ overlap and inter-W BEC, a 
detailed analysis of BEC in a single W would be required. Given that only 
a few thousand of these W's have been produced at LEP, such a detailed 
study is presently not possible. It is, however, possible on the millions of 
events accumulated at the Z, and we shall assume that the fragmentation 
properties are similar for those two bosons (except for the fact that Z — > bb 
has no equivalent in W decay) . Since this detailed analysis is still going on, 
we shall also look at corresponding properties in hadron-hadron and even 
heavy-ion collisions. 

4-2. Pion-source elongation 

The form of the correlation function in more than one dimension has 
been a major subject of theoretical study in recent years[28, 34, 35, 19, 
36, 37, 38]. In Monte Carlo generators, spherical symmetry is usually 
assumed[5, 39, 40, 41], while elongation can be expected when a string- 
like shape is maintained [28, 38]. Experimentally, detailed three-dimensional 
analyses were done for heavy- ion collisions [42, 43] and for hadron-hadron 
collisions [44]. While the volume of the pion emission region appeared to 
be approximately spherical for heavy-ion collisions, a clear elongation was 
observed in hadron-hadron collisions. An elongation is now also observed 
at LEP [45, 46, 47]. 
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Figure 4. The longitudinal center of mass frame (LCMS) showing the projection of 
Q on the (QL-Qout) plane. Q S idc is the projection of Q on the axis perpendicular 
to this plane. 



In this analysis the longitudinal center-of-mass system (LCMS) [36] is 
used. This is defined for each pair of particles as the system, resulting from 
a boost along the thrust axis, in which the sum of the momenta of the 
pair is perpendicular to the thrust axis. In this system, one can resolve the 
three-momentum difference of the pair of particles into a longitudinal com- 
ponent Ql parallel to the thrust axis, Qout along the sum of the particles' 
momenta (see Fig. 4) and Q s ide perpendicular to both Ql and Qout- Then, 
the invariant four-momentum difference can be written as [36] 



Q 2 = Ql + Qside + Qout " (A£) 2 = Ql + QLe + Qout(l " A (27) 

where 

o _ Pout 1 ~\~ Pont 2 foe"! 

E 1 + E 2 1 ' 

with p ou ti and (i = 1,2) the out-component of the momentum and the 
energy, respectively, of particle i in the LCMS. The energy difference AE 
and therefore the difference in emission time of the two particles couples 
only to the component Qout- Consequently, Ql and Q s ide reflect only spa- 
tial dimensions of the source, whereas Qout reflects a mixture of spatial 
and temporal dimensions. The correlation function is then parametrized in 

terms of Q = (Ql, Qside, Qout): 



(29) 
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Assuming a Gaussian (azimuthally, but not necessarily spherically, sym- 
metric) shape of the source, the following three-dimensional parametrization 
has been proposed [35, 19, 48]: 

RliQhi Qoutj Qsidc) = 

= 7 (1 + SQl + eQout + £Qsidc) • (30) 
• [l+Aexp (^rlQl-rl nt Ql ut -rl idc Ql iAc +2p^ out r^r out Q h Q out )] , 

where the factor (l + 5QL + £<3out+?Qside) takes into account possible long- 
range momentum correlations in the form of a slow rise, 7 is a normalization 
factor close to unity and the term between square brackets is the two-particle 
Bose-Einstein correlation function associated with a Gaussian shape of the 
source. 

By fitting the correlation function with this parametrization, one can ex- 
tract the factor A, which measures the strength of the correlation, and the 
'radii' vi (i = L, out and side) defined as ai/y/2, with the of the variances 
of a multi-dimensional Gaussian distribution of the source in configuration 
space. pL,out is the correlation between the longitudinal and out components 
of this Gaussian. In the LCMS, the duration of particle emission only cou- 
ples to the out-direction and only enters in the parameters r out and PL,out- 
Hence, r s id e can be interpreted as the transverse component of the geomet- 
ric radius. The parametrization, Eq. (30), assumes azimuthal symmetry 
of the source, which means that the two-particle Bose-Einstein correlation 
function associated with the Gaussian shape of the source, is invariant un- 
der the transformation Q s ide — ► — Qside- Consequently, the only possible 
off-diagonal term is the QhQout term. It turns out to be zero within errors, 
however. Note that /9L,out = is indeed expected in LCMS near midrapidity 
or for boost-invariant sources. 

The results of three LEP experiments [45, 46, 47] are summarized in 
Table 1. In spite of the different selection criteria and reference samples, all 
experiments consistently demonstrate an elongated shape of the pion source 
(or rather region of homogeneity) in hadronic Z decay. On the other hand, 
r s idc/r L = 1-08 ± 0.03 is found [45] for JETSET with BE [49]. 

A systematic study of the hierarchy of radii obtained from JETSET was 
recently performed in [50] . Starting from a spherically symmetric Gaussian 
correlator, the authors obtain 

^side > r L > r out , 

unlike the experimentally observed elongation, both when using momentum 
shifting [49] or event weighting [51] to simulate the correlation. Generaliz- 
ing to asymmetric weights, the experimentally observed elongation can be 
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Table 1. Elongation of the pion source in hadronic Z° decays: summary of the 
measurements at LEP1 (rx corresponds to Qt = yQZut + Qsidc)- 





L3 

"mixed" reference, 
all events 


DELPHI 

"mixed" reference, 
2- jet events 


OPAL 

"H — " reference, 
2-jet events 


A 


o.4i±o.oilS:SS 


0.261 ±0.007 ±0.010 


0.443 ± 0.005 


r-L, fm 


0.74 ± 0.021°;°^ 


0.85 ±0.02 ±0.07 


0.989 ±0.0111q.o?b 


rout, fm 


0.53 ± 0.021°;°^ 




0.647 ±0.011^24 


r s idc, fm 


0.59 ± O.OllHs 




0.809 ± 0.009io;°32 


r s ide/r L 


0.80 ± 0.02±g;?| 




0.818 ± 0.018±g;ggS 


tt, fm 




0.53 ±0.02 ±0.07 




r T /r L 




0.62 ±0.02 ±0.05 





reproduced, but finding a good set of input radius parameters turns out an 
involved procedure. 

What is important to realize is that the measured longitudinal radius 
has nothing to do with the elongation of the qq string stretched in Z decay. 
We shall see in Sect. 4.6 that the full pion-emission function is of the order 
of 100 fm long, so that tl only measures a small fraction of it. The reason 
for that is a strong x, p correlation. Pions produced at a large distance 
on the string also have very different momenta and do not correlate. The 
"radii" therefore measure the effective size of the source segment radiating 
mesons with sufficiently small relative momentum (length of homogeneity) , 
as shown in Fig. 5. 




z 



Figure 5. Space-time picture of particle emission for a given fixed mean momentum 
of the pair. The mean value of the proper-time and the space-time rapidity distri- 
butions is denoted by f and rj. As the rapidity of the produced particles changes 
from the smallest to the largest possible value, the \f(y),r](y)] variables scan the 
surface of mean particle production in the (t, z) plane [38] . 
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4-3. The functional form of the correlation function 



More important than the parameters extracted from "forcing" the two- 
particle correlation function into a fit by a pre-selected parametrization, is 
the actual experimentally observed shape of this distribution, itself. 

The simple geometrical interpretation of the interference pattern based 
on the optical analogy as in Sect. 3.1 is invalid when emitters move relativis- 
tically with respect to each other, leading to strong correlations between the 
space-time and momentum-energy coordinates of emitted particles [52, 53]. 
Correlations of this type arise due to the nature of inside-outside cascade 
dynamics [54] as in colour-string fragmentation [55]. In the interpretation of 
BEC by Andersson and Hofmann [23] in the string model, the length scale 
measured by BEC is therefore not related to the size of the total pion emit- 
ting source, but to the space-time separation between production points 
for which the momentum distributions still overlap. This distance is, in 
turn, related to the string tension. The model predicts an approximately 
exponential shape of the correlation function 



where r is expected to be independent of the total interaction energy. 

Furthermore, scale-invariant dynamics is strongly connected with Bose- 
Einstein correlation. Scale invariance implies that multiparticle correlation 
functions exhibit power-law behavior over a considerable range of the rel- 
evant relative distance measure (such as Q 2 ) in phase space [6, 7]. As 
such, BEC from a static source do not exhibit power-law behavior. How- 
ever, a power law is obtained if the size of the particle source fluctuates 
event-by-event, and/or, if the source itself is a self-similar (fractal-like) ob- 
ject extending over a large volume [56]. In these studies, the ratio R2 is 
parametrised using the form 



The usually "reasonable" \ 2 values of the Gaussian fits hide the fact 
the Gaussian parametrization in general fails at low values of Q 2 , where 
statistical errors are often large. For the case of two-particle correlations, 
this has been demonstrated convincingly by NA22 [57] and UAl [58], but 
deviations from a Gaussian are also observed in lepton-hadron [59, 60] and 
e + e~ [30] collisions. 

This failure of a (multi-) Gaussian form persists in higher-order correla- 
tions. In Fig. 6a the NA22 data [61] on BE correlations of order q = 2 to 4 
are plotted as a function of Q 2 , in conventional linear scale. The curves are 



R 2 (Q) = R {1 + Aexp(-rQ)) , 



(31) 




(32) 
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Figure 6. The normalized two-, three- and four-particle inclusive densities as a 
function of Q 2 (left) and -InQ 2 (right) [61]. Curves show the multi-Gaussian fits 
according to [62]. 



the fits by a (/-fold Gaussian parametrization [62]. In Fig. 6b the same data 
and the same fits are repeated for Q 2 < 1 GeV 2 on ln-ln scale (where the Q 2 
axis is reflected, i.e., small Q 2 correspond to large — InQ 2 ). Even though the 
statistical errors at small Q 2 are large (the very reason why small Q 2 does 
not contribute much to x 2 ), it is obvious that small-Q 2 points systematically 
lie higher than the multi-Gaussian fit, thus supporting a power-law behav- 
ior. This effect is even enhanced when the data are corrected for Coulomb 
repulsion. 

Fig. 7a shows [63] the second-order cumulant if 2 as a function of Q 
(on log-scale) compared to a general quantum statistical model, based on 
a classical source current formalism applied successfully in quantum optics 
[64]. It includes as special cases more specific models such as [18] and 
[62] . The APW normalized cumulant predictions are built from normalized 
correlators dij , the on-shell Fourier transforms of classical space-time current 
correlators. The specific parametrizations tested in Fig. 7 are 

Gaussian : = exp(— r 2 Q 2 ) 
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Q, GeV max(Q 12 ,Q 2 3,Q 3 i), GeV 

Figure 7. a) Second-order cumulant with fits of the forms given, b) Third-order 
cumulant with APW predictions based on K 2 [63]. 



exponential : dij = exp(— rQij) (33) 
power law : = Q~ Q . 

For constant chaoticity A and real- valued currents, APW predict as second- 
and third-order normalized cumulants 

K 2 {Q 12 ) = 2A(1 - A)dia + \ 2 d\ 2 , (34) 

K 3 (Qi2, Q23, Q31) = 2A 2 (1 - X)[d 12 d 23 + d 23 d 31 + d 31 d 12 ] + 2X 3 d 12 d 23 d 31 . 

(35) 

The fits in Fig. 7a contain an additive "background" parameter, in addition 
to A and r as free parameters. The Gaussian fit is clearly excluded and the 
best fit is achieved with the power-law parametrization of the correlation 
function (full line). 

Also K 3 plotted in Fig. 7b shows a power-law increase. It is, furthermore, 
visible in Fig. 7b that the increase is faster than expected from APW K 2 , 
even for the power-law parametrization. 

So, there is ample room for improvement of the models and we believe 
that the recently developed methods of studying the correlations (higher- 
order cumulants, higher dimensionality, alternative parametrizations of the 
correlation function) have opened the way for an improvement of these 
models. 

An interesting extension of the usual Gaussian approximation of the BE 
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Figure 8. Projection onto Qt of a two-dimensional Gaussian (dashed) and Edge- 
worth (solid) fits for a small-QL slice [65] . 

correlation function is an Edgeworth expansion [16] as suggested in [65] 
R 2 (Q) = 7 (1 + A* exp(-t 2 /2)[l + ^H 3 (t) + ^H 4 (t) + ...]), (36) 

with t = V2Q • r, H n being the n-th Hermite polynomial, and K n the n-th 
order cumulant moment of the correlation function, where k<i yields r. The 
Hermite polynomials of odd order vanish at the origin, so that 

A = A*[l + /c 4 /8 + ...] . (37) 

A generalization to higher dimensions is straightforward [65], except for 
possible correlations between the Qi variables. 

The influence of the non-Gaussian shapes was studied [65] on AFS [66, 
67], E802 [68] and NA44 [69] data. In Fig. 8, a Q T projection of a 2D 
Edgeworth fit is compared to that of a 2D Gaussian fit to the E802 data. 
The deviation from a Gaussian (dashed) is obvious, and the Edgeworth 
expansion (full line) is flexible enough to describe it (with A = 1!). 

In Fig. 7a it was shown that even an exponential is not steep enough 
to reproduce the fast increase of K 2 - An interesting observation of [65] is, 
that a Laguerre expansion of an exponential can reproduce these UAl and 
the NA22 data (Fig. 9). However, at low Q 2 data are still systematically 
above the fit and a power-law fit is reported in [65] to give similarly good 
X 2 /NDF with a smaller number of fit parameters. With a core-halo model 
[70] strength parameter of A* = 1.14 ± 0.10 (UAl) and A* = 1.11 ± 0.17 
(NA22), i.e., at maximum possible value (unity), there are either other than 
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BE correlations at work or all resonances are resolved at these low Q 2 values. 
This may imply the connection between the observed power-law behavior 
(intermittency) and resonance contributions of BE correlations [56] . 



Lo guerre expansion fit 




Figure 9. The figures show F| which is proportional to the two-particle Bose- 
Einstcin correlation function, as measured by the UA1 [58] and the NA22 [57] 
Collaborations. The dashed lines stand for the exponential fit, the solid lines for 
that with the Laguerre expansion [65] . 

In a 3D analysis of e + e~ collisions at the Z-mass [45], results more satis- 
factory than those obtained with either Gaussian or exponential parametriza- 
tions were obtained with the Edgeworth expansion. Taking only the lowest- 
order non-Gaussian term into account (and dropping the off-diagonal term) , 
Eq.(30) results in 

R2(Qh, Qoutj Qsidc) = 

7(1 + 5Q L + eQ out + CQside) 



■{l + Xexp(-rlQt 



3! 

where Ki (i 



1- 



out *e out 
K out 



r sidcQside) 



(38) 

H 3 (r utQont) l+^|p^3(^sideQside)] } , 

L, out, side) is the third-order cumulant moment in the corre- 
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0.65 ± 0.02 


0.91 ±0.01 


OPAL [77] 
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1.06 ±0.05 ±0.16 


0.49 ±0.01 ±0.05 


DELPHI [30] 


prompt 
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K±K± 


0.82 ±0.11 ±0.25 


0.48 ± 0.04 ± 0.07 


DELPHI [75] 


all 


unlike 




0.82 ± 0.22l°;^ 


0.56 ± 0.08t° -° l 


OPAL [76] 


2-jet 


mixed 


KgKg 


1.14 ±0.23 ±0.32 


0.76 ±0.10 ±0.11 


OPAL [78] 


all 


MC 




0.61 ±0.16 ±0.16 


0.55 ±0.08 ±0.12 


DELPHI [75] 


all 


MC 




0.96 ±0.21 ±0.40 


0.65 ±0.07 ±0.15 


ALEPH [79] 


all 


MC 



Table 2. Parameters Aq and tq in the Gaussian parametrization in e + e interac- 
tions at LEP, for different like-charged particles [76] . 

sponding direction and H^(RiQi) = (V^RiQi) 3 —Sy/2RiQi is the third-order 
Hermite polynomial. Note that the second-order cumulant corresponds to 
the radius r^. Applying this expansion to the L3 data [45] discussed in 
Sect. 4.2 improves the confidence level of the fit from 3% to 30%. Non- 
zero values of the k parameters indicate the deviation from a Gaussian, A is 
larger than the corresponding Gaussian A and the values of the radii confirm 
the elongation observed from the Gaussian fit. 

4-4- (Transverse) mass dependence 
44.1. The K^K^ system 

Kaons are less affected by resonance decay than pions and could eventu- 
ally provide a cleaner signal of the source. Bose-Einstein correlations among 
equally-charged kaons were observed in hh [71, 72], AA [73, 74] and e + e~ 
[75, 76] collisions (see Table 2 for the latter). 

The size of the kaon emission region tends to be smaller than that of 
the pion emission region, in particular in AA collisions. The difference 
in resonance effects on tttt and KK correlations only partially can explain 
this difference in AA collisions. In e + e~ collisions, the Gaussian radius 
parameter re tends to be smaller for K ± K ± than for it^it^, but the spread 
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is large due to different choices of background. 

44.2. The K° S K° S system 

The K g K s system is a mixture of K°K° and K°K° (K°K°) pairs. At 
LEPI energy, only 28% of all KgKg pairs are estimated to come from the 
(identical) K°K° or K°K° system. What is particularly interesting is that 
Kg's can interfere even if they originate from a (non- identical) K°K° system 
[80]: An enhancement is expected in the low-Q region if one selects the 
C = +1 eigenstate of 

|K°K° > c=±1 = i=(|K°(p)K°(-p)> ± K°(p)K°(-p))), (39) 

where p is the three-momentum of one of the kaons in their cms. In the limit 
p -► (Q -► 0), the C = -1 (KgK£) state disappears and C = +1 (KgK s 
or K^Kl) becomes maximal. 

The enhancement in KgKg and K^K^ pairs at low Q is exactly compen- 
sated by the low Q suppression of the KgK^ state, so that no BE effect is 
to be expected as long as all possible final states of the K°K° system are 
considered. A full BE-like enhancement is, however, expected for the KgKg 
system by itself. 

Early, low statistics results come from the hh experiment [81], new re- 
sults exist from DELPHI [75, 82], OPAL [78] and ALEPH [79]. While the 
kaon-production radius is smaller for the hh experiment, it seems to agree 
with those measured for both charged kaons and pions in the e + e~ ex- 
periments, within the large spread of values observed. Furthermore, the 
parameter A is large in agreement with the expectation [80]. 

44.3. A A or A A 

An interesting generalization of the Bose-Einstein formalism used above 
is to consider Fermi-Dirac interference, essentially by changing the sign in 
front of the correlator. This leads to a destructive interference at small 
phase-space distance and allows to determine the emission radius for iden- 
tical fermions in a comparison of the amount of their total-spin 5 = 1 state 
(destructive) to that of their S = state (constructive) as a function of Q 
[83]. The method does not need a further reference sample. It was applied 
to e + e~ data at LEPI in [84, 85, 86] and gives a radius of about 0.15 fm. 
It was, however, verified, that the conventional method with JETSET as a 
reference sample gives similarly low a radius. 
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4-4-4- (Transverse) mass dependence of the radius parameter 

The simultaneous comparison of the emission radii for pions, kaons and 
A's now suggests a decrease with increasing mass. Such a behavior has 
first been observed by NA44 in heavy- ion collisions [73]. The NA44 results 
can be translated into an l/^/m T scaling of the radius, in agreement with 
the expectations from a hydrodynamical model [87] with three-dimensional 
collective expansion and cylindrical symmetry. 

In Fig. 10a, the radius parameter r is shown as a function of the hadron 
mass m [88] for e + e~ annihilation at the Z mass. The large error associated 
with r n7r reflects the systematic uncertainty due to the choice of the reference 
sample. A general trend can be observed as a hierarchy 

TVtt > r KK > r AA ■ (40) 



Some effect is to be expected from kinematics, i.e. from the mass-dependent 
integration limits when transforming from i?2(Pi>P2) i n six-dimensional 
momentum space to B,2(Q) in one-dimensional momentum separation [89]. 
This effect is far too small, however. 

The authors [88] show that a \jyfm behavior can be expected already 
from the Heisenberg principle with 

ApAr = mvr = he 

AEAt = p 2 At/m = h (41) 
cVhAt 

and r = — 7 =^- , (42) 
sjm 

where m, v and p are the hadron mass, velocity and momentum and r is the 
distance between the two hadrons. Assuming AE to only depend on the 
kinetic energy of the produced particle and At = 10~ 24 sec, independent of 
m, grants the thin solid line in Fig. 10a. The upper and lower dashed lines 
correspond to an increase or decrease of At by 0, 5 • 10~ 24 sec, respectively. 
(The thick solid line corresponds to a perturbative QCD cascade using the 
virial theorem and assuming local parton hadron duality (LPHD).) 

However, as shown in [88], a formula identical to (42) also holds for 
the radius r z in the longitudinal direction and the average transverse mass 
m T = 0.5 (y/m 2 +p2, i + ^J m 2 +p ^ 2 ). Fig. 10b shows DELPHI results [90] 

compared to At = 10~ 24 sec (dashed) and the best fitted value of At = 
2.1 • 10~ 24 sec (full line). 

Alternatively, the transverse mass dependence can be explained by a 
generalized inside-outside cascade [91] assuming (i) approximate proportion- 
ality of four-momenta and production space-time position (freeze-out point) 
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Figure 10. a) The radius parameter rasa function of the hadron mass m. b) The 
longitudinal emitter radius r z as a function of mi [90]. The lines are described in 
the text [88] . c) Longitudinal and transverse radius as a function of the transverse 
mass M± of the two-particle system, compared to the Mj_-threshold values (same 
as data in Fig. 10a) [91]. 



of the emitted particles = ax^ and (ii) a freeze-out time distributed along 
the hyperbola Tq = t 2 — z 2 (i.e., a generalization of the so-called Bjorken- 
Gottfried conditions). Prom the two conditions above follows directly 

a 2 r 2 = E 2 - p 2 z = m\ (43) 

and the generalized Bj or ken- Gottfried condition 



772-T 

Pu = x u 

TO 



(44) 
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Using a more rigorous formulation in terms of the Wigner representation, 
the authors show how this proportionality leads to an wit dependence of 
the radius parameter. Fig. 10c gives indeed a dependence of both the longi- 
tudinal and the transverse radius on the transverse mass of the two-particle 
system, 

Ml=fh\ + m T1 m T 2 sinh 2 ^ ~ 2/2 ^ . (45) 

For a set of "reasonable" model parameters [91], the experimental results 
(same as in Fig. 10a) are reproduced reasonably well. Note that the exper- 
imental data are given at the threshold value of the corresponding M_l at 
which transverse momenta and rapidity differences are small compared to 
the particle masses. 

The parameters are to be improved, but is closely related to the av- 
erage transverse momentum and S± has to be considerably larger than 
in the model to satisfy the uncertainty principle. Since S± corresponds to 
a correlation length between transverse momentum and transverse position 
at freeze-out, this correlation is rather weak. Nevertheless, it is sufficient to 
create a strong variation of the transverse radius, and suggests the existence 
of an important "collective flow" , even in the system of particles produced 
in e + e~ annihilation! Note that strong space-time momentum-space corre- 
lations are expected not only from hydrodynamic expansion, but also from 
jet fragmentation. 



4-5. The multiplicity (or density) dependence 

In nucleus-nucleus experiments [92, 93], the radius r was found to in- 
crease with increasing charged-particle multiplicity n. By relating r to the 
size of the overlap region of the two colliding particles, this increase can be 
understood in terms of the geometrical model [94] : a large overlap should 
imply a large multiplicity. On the other hand, no evidence for a multiplicity 
dependence is found in hadron-nucleus collisions at 200 GeV/c [95]. 

After some time of confusion, the n dependence is now clear for hadron- 
hadron collisions. At energies below y/s pa 30 GeV (i.e at y/s pa 8 [96], 
22 [97] and 27 GeV [72]) no n-dependence is observed for re- At higher 
energies (last ref.[98] and [71]) an n-dependence starts to set in and to grow 
with increasing energy (see Fig. 11a). At the highest ISR energy (y/s = 62 
GeV) the increase is about 40% when the density in rapidity is doubled, 
but at \/s = 31 GeV the increase is still very weak. The result is extended 
to yfe = 630 GeV by UA1 [99] and to 1800 GeV by E735 [100] in Fig. lib. 
At very large density, the increase of tq with increasing density is shown 
to extrapolate well to the heavy- ion results of NA35 [93] in Fig. 11c. The 
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Figure 11a) Radius tq of the pion source as a function of charged-particle density 
for the energies indicated [98], b) same for pp collisions at 630 GeV [99] and 1800 
GeV [100], as well as e + p collisions at 300 GeV [60], c) Comparison of tq as 
a function of charged-particle density An/Arj [99] with the results of relativistic 
heavy- ion collisions [93]. 



effect is reproduced in thermodynamical and hydro dynamical models. The 
A parameter, on the other hand, decreases with increasing n (not shown). 

At the low-density side, the effect is also observed in e + p collisions by 
HI [60] (crosses in Fig. lib). The results from e + e~ experiments at lower 
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Figure 12. a) Dependence of Aq and re on the charged-particle multiplicity n for 
c+c~ collisions at the Z mass, b) same for two-jet events (solid points) and three-jet 
events (open points) [77]. 

energy [26, 101] were consistent with no multiplicity dependence as expected 
from the geometrical model, but also for this type of collisions a multiplicity 
dependence was finally established at higher energy [77] (see Fig. 12a). At 91 
GeV, the radius vq is found to increase linearly with increasing multiplicity 
n, showing a small but statistically significant increase of about 10% for 
10 < n < 40. As for hh-collisions, the chaoticity parameter Ag decreases 
with increasing n. 

In Fig. 12b, OPAL further shows that the multiplicity dependence is 
strongly reduced in separate samples of two-jet and three-jet events, the 
average value of vq, however, being 10% bigger for three-jet than for two-jet 
events. Folding in the multiplicity difference of two- and three-jet events, 
this at least partly explains the effect as due to multi-jet production at 
higher energies. The decrease of A is larger in the 3-jet than in the 2-jet 
sample. 

As shown quantitatively in [102], it is crucial to study the normalized 
cumulants ^(Q) (Eq.(ll)) rather than the normalized densities R-ziQ) 
(Eq.(10)) in a density dependent analysis and to correct for a well-defined 
multiplicity-dependent bias due to the cut in the multiplicity distribution 
(the point being that K q / for limited n, even in case of independent 
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emission). In Fig. 13a) and b) [103], the bias-corrected (so-called "inter- 
nal") cumulants are given for UA1 as a function of the inverse rapidity 
density, for small and large values of Q, respectively. The data show 

i) a linear dependence (similar for like-charged and unlike-charged pairs), 

ii) vanishing of the cumulant at large density for large Q, 

iii) approach towards a finite limit for large density at small Q (where 
BE correlations are expected to dominate). 




1/(dN c /d7j) 1/(dN c /d7j) 



Figure 13. Inverse density dependence of the bias-corrected cumulant at Q = 0.1 
GeV (Fig. 13a) and Q = 7.0 GcV (Fig. 13b) for like-charged pairs (full circles) and 
unlike-charged pairs (open circles). The crosses in Fig. 13a correspond to A- values 
[103]. 

The large-Q behavior [points i) and ii) above] is that expected from 
particle emission from N fully overlapping, identical but fully independent 
sources (e.g. strings). From the additivity of unnormalized cumulants fol- 
lows immediately a dilution, 

^ ■ ■ ■ •*> = - ■■■■*>■ ^ 

where Kq N ^ is the normalized cumulant of the iV-source system, while K^p 
is that of an individual source (see also [104, 105, 106]). This results in 
a normalized cumulant inversely proportional to N or to the total density 
dn/dy = Ndn^ /dy, as observed in Fig. 13b. 

At small Q (Fig. 13a), however, the normalized cumulant approaches a 
constant different from zero at large densities. Naively, this would imply 
correlations between particles coming from different sources, which could be 
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interpreted as inter-source Bose-Einstein correlations, would not a similar 
effect be observed for unlike-charged pairs, as well. So also resonances play 
an important part. One has to keep in mind, however, that (46) only holds 
for full overlap of identical sources and that at Q ~ 0.1 GeV the overlap is 
far from complete and the number of sources limited. 

However, also in heavy-ion collisions there is evidence that A does not 
drop with increasing density for ever. Heavy-ion collisions lead to A values 
quickly decreasing with increasing density at lower densities (i.e. lower- A 
collisions). In agreement with the expectation for overlapping independent 
sources, A drops from 0.79 to 0.32 [107] from O-C to O-Cu, O-Ag and O- 
Au. In high- A collisions, on the other hand, a saturation seems to set in 
[108, 109]. For S-Pb and Pb-Pb central collisions NA44 [73] quotes A = 0.56 
and 0.59, respectively. Such a saturation and eventual increase of A would 
be expected if the densely packed strings of a heavy-ion collision finally 
coalesce until they form a large single fireball (percolation of strings) [109]. 



As has become clear from the previous sub-sections, the correlation mea- 
surements alone do not contain the complete information on the geometrical 
and dynamical parameters characterizing the evolution of the hadronic mat- 
ter. In particular, BEC are not measuring the full geometrical size of large 
and expanding systems, since that expansion results in strong correlations 
between space-time and momentum space. More comprehensive information 
can be provided by a combined analysis of data on two-particle correlations 
and single-particle inclusive spectra [70, 21, 110, 111, 112]. 



In the framework of the hydrodynamical model for three-dimensionally 
expanding cylindrically-symmetric systems [70] , the emission function corre- 
sponds to a Boltzmann approximation of the local momentum distribution. 
Within this model, the invariant single-particle spectrum of pions in rapidity 
y and transverse mass mx is approximated by 



4-6. The emission function 



4-6.1. The formalism: 



f{y,m T ) = 



1 dN w 



C„ !T Vo S h, fceX p(M)exp[-<^£ 




( m T \ 
exp^--jx 




2T [T + ((u T y + (^-))m T ] 



(u T } 2 (m^ - ml) 



(47) 
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with 



Ay* = A V 2 + ^- (48) 

7BT 

1 1 TUT 

+ ~TrT~ coshr/ s , (49) 



At/ 2 At/ 2 T 

y - yo 



1 + At/ 2 ^" 



(50) 



The width Ay of the rapidity distribution given by (48) is determined by 
the width At/ of the longitudinal space-time rapidi ty t/ dis tribution of the 
pion emitters and by the thermal smearing width y/T^J ttit , where To is the 
freeze-out temperature (at the mean freeze-out time Tf) at the axis of the 
hydrodynamical tube, To = Tf(rx = 0). For the case of a slowly expanding 
system one expects At/ <C To/mi, while for the case of a relativistic lon- 
gitudinal expansion the geometrical extension At/ can be much larger than 
the thermal smearing (provided tot > To). 

In addition to the inhomogeneity caused by the longitudinal expansion, 
(47) also considers the inhomogeneity related to the transverse expansion 
(with the mean radial component (wt) of hydrodynamical four- velocity) and 
to the transverse temperature inhomogeneity, characterized by the quantity 



/ AT\ JTq_ 
\ T / T r 



- 1, (51) 



mis 



where T rms = Tf(rr = rT(rms)) is the freeze-out temperature at the trans- 
verse rms radius rT(rms) and at time Tf. 

The power a in (47) is related [70] to the number d of dimensions in 
which the expanding system is inhomogeneous. For the special case of the 
one-dimensional inhomogeneity (d = 1) caused by the longitudinal expan- 
sion, a = 1 — 0.5d = 0.5 (provided At/ 2 3> To/tot). The transverse inho- 
mogeneity of the system leads to smaller values of a. The minimum value 
of a = —1 is achieved at d = 4 for the special case of a three-dimensionally 
expanding system with temporal change of local temperature during the 
particle emission process. 

The parameter uq in (47) denotes the midrapidity in the interaction 
c.m.s. and can slightly differ from due to different species of colliding 
particles. The parameter C is an overall normalization coefficient. 

Note that (47) yields the single-particle spectra of the core (the cen- 
tral part of the interaction that supposedly undergoes collective expansion). 
However, also long-lived resonances contribute to the single-particle spectra 
through their decay products. Their contribution can be determined in the 
core-halo picture [87, 113] by the momentum dependence of the strength 
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parameter A(y, hit) of the two-particle Bose-Einstein correlation function. 
Experimentally, the parameter is however found to be approximately inde- 
pendent of tot [114, 115, 116]. Hence this correction can be absorbed in 
the overall normalization. 

The two-dimensional distribution (47) can be simplified for one-dimen- 
sional slices [70, 38, 116]: 

1. At fixed mT, the rapidity distribution reduces to the approximate 
parametrization 



f(y,m T ) = C m exp 



2Ay 2 



(52) 



where C m is an mi-dependent normalization coefficient and yo is defined 
above. The width parameter Ay 2 extracted for different mx-slices is pre- 
dicted to depend linearly on l/m-p, with slope To and intercept Arf (cf. 
(48)). 

Note, that for static fireballs or spherically expanding shells (52) and 
(48) are satisfied with Arj = [38]. Hence, the experimental determination 
of the 1 / tot dependence of the Ay parameter can be utilized to distinguish 
between longitudinally expanding finite systems versus static fireballs or 
spherically expanding shells. 

2. At fixed y, the m T -distribution reduces to the approximate parametri- 
zation 

f(y,m T ) = C y m%exp (~^) > ( 53 ) 

where C y is a y-dependent normalization coefficient and a is defined as 
above. 

The y-dependent "effective temperature" T e g(y) can be approximated 

as 

I± 

1 + a(y - y f 



T M = ; , , 2 , (54) 



where T* is the maximum of T e g(y) achieved at y = yo, and 

a = T ° T * (55) 

2m 2 (Ar?2 + JL )2 

with To and Arj 2 as defined above. 

The approximations (52) and (53) explicitly predict a specific narrowing 
of the rapidity and transverse mass spectra with increasing tot and y, re- 
spectively (cf. (48) and (54)). The character of these variations is expected 
[38] to be different for the various scenarios of hadron matter evolution. 
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4-6.2. The results: 

The Ay 2 values obtained from fits of the NA22 data [116] by (52) are 
given as a function of 1/rriT in Fig. 14a. A fit to the widening of the rapidity 
distribution (i.e. increase of Ay 2 ) with increasing 1/mx by (48) gives an 
intercept Arj 2 = 1.91±0.12 and slope T = 159±38 MeV. Thus, the width of 
the y-distribution is dominated by the spatial (longitudinal) distribution of 
pion emitters (inherent to longitudinally expanding systems) and not by the 
thermal properties of the hadron matter, as would be expected for static or 
radially expanding sources. Since Arj 2 is significantly bigger than 0, static 
fireballs or spherically expanding shells, able to describe the two-particle 
correlation data in [115], fail to reproduce the single-particle spectra. 
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Figure 14. a) The (l/mT)-dependcncc of (Ay) 2 for inclusive ir~ meson rapidity dis- 
tributions at | y | < 1.5. The straight line is the fit result according to paramctriza- 
tion (48). b) T e g as a function of y fitted according to parametrization (54) [116]. 



The T c g values obtained from fits of the same data by (53) are given as 
a function of y in Fig. 14b. T e g (y) tends to decrease with increasing \y\ and 
approximately follows (54) with T* = 160 ± 1 MeV, a = 0.083 ± 0.007 and 
yo = —0.065 ± 0.039. Note, however, an asymmetry in the T e Q distribution 
with respect to y = 0: except for the last point, T e g is higher in the meson 
than in the proton hemisphere. 

The values of the exponential parameter a fitted in (53) are near zero, 
corresponding to a two-dimensional inhomogeneity of the expanding system 
(a = 1 — 0.5<i) . One concludes, therefore, that apart from a longitudinal in- 
homogeneity caused by the relativistic longitudinal flow, the hadron matter 
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also has a transverse inhomogeneity (caused by transverse expansion or a 
transverse temperature gradient) or undergoes a temporal change of local 
temperature during the particle emission process. 



4-6.3. The transverse direction: 

Further information on hadron-matter evolution in the transverse direc- 
tion can be extracted from (47) with parameters (ut) and (■ A f-) characteriz- 
ing the strength of the transverse expansion and temperature inhomogene- 
ity. 

A moderate value of the mean transverse four-velocity (ut) = 0.20 ± 
0.07 indicates that the transverse inhomogeneity is mainly stipulated by the 
rather large temperature inhomogeneity {^-) = 0.71 ± 0.14. Using (51), 
one infers that the freeze-out temperature decreases from Tq = 140 ±3 MeV 
at the central axis of the hydrodynamical tube to T rms = 82 ± 7 MeV at a 
radial distance equal to the transverse rms radius of the tube. 



4.6.4. Combination with two-particle correlations: 

Due to the non-static nature of the source, the effective size parameters 
r L ; '"out j '"side vary with the average transverse mass wt = ^(^ti + m^) 
and the average rapidity Y = \{yi + 2/2) of the pion pair. In the LCMS the 
effective radii can be approximated [70, 38, 111] by 

r 2 h = r 2 A V 2 . (56) 

r 2 ont =r 2 ,+fcA T 2 (57) 



with 



lite = < (58) 



rl 



r * 1 1 mi// \2 1 /AT\\ ' 

where the parameters An 2 , To, («t) an d {^-) are defined and estimated 
from the invariant spectra above; r g is related to the transverse geometrical 
rms radius of the source as r s (rms) = V%r g ; Tf is the mean freeze-out 
(hadronization) time; Ar* is related to the duration At/ of pion emission 
and to the temporal inhomogeneity of the local temperature. If the latter 
has a small strength (as one can deduce from the restricted inhomogeneity 
dimension estimated above), an approximate relation At/ > At* holds. 
The variable /?t is the transverse velocity of the pion pair. 
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Far am. 


NA22 


TT T 

Heavy ion 
Averaged 


T [MeV] 


140 ±3 


139 ±6 


(u T ) 


0.20 ± 0.07 


0.55 ± 0.06 


r g [fm] 


1.2 ±0.2 


7.1 ±0.2 


Tf [fm/c] 


1.4 ±0.1 


5.9 ±0.6 


At/ [fm/c] 


1.3 ±0.3 


1.6 ± 1.5 


At? 


1.36 ± 0.02 


2.1 ±0.4 


<^> 


0.71 ±0.14 


0.06 ± 0.05 


yo 


0.082 ± 0.006 


(fixed) 



Table 3. Fit parameters of the Buda-Lund hydro (BL-H) model in a combined 
analysis of NA22 [116], NA49, NA44, WA98 [117] spectra and correlation data. 

Using (56) and (59) with T = 140 ± 3 MeV and At? 2 = 1.85 ± 0.04, 
together with tl fitted in different % ranges, one finds a mean freeze-out 
time of Tf = 1.4 ± 0.1 fm/c. 

The transverse-plane radii r out and r s id c measured in [115] for the whole 
ttt-t range are: r out = 0.91 ±0.08 fm and r s jd e = 0.54±0.07 fm. Substituting 
into (57) and (58), one obtains (at fh = 0.484c [115]): At* = 1.3 ± 0.3 
fm/c. Since the mean duration time of pion emission can be estimated as 
At/ > At*, the data grant At/ Tf. A possible interpretation is that in 
meson-proton collisions the radiation process occurs during almost all the 
hydrodynamical evolution of the hadronic matter produced. 

An estimation for the parameter r g can be obtained from (58) and (60) 
using the quoted values of r s id e ,?o, (ut) and {^-). The geometrical rms 
transverse radius of the hydrodynamical tube, r 9 (rms) = V2r g = 1.2 ± 0.2 
fm, turns out to be larger than the proton rms transverse radius. 

The set of parameters of the combined analysis of single-particle spectra 
and Bose-Einstein correlations in (-7r + /K + )p collisions [116] is compared to 
that obtained [117] from averaging over Pb±Pb experiments (NA49, NA44 
and WA98) in Table 3. 

The temperature To near 140 MeV comes out surprisingly similar for hh 
and PbPb collisions. The geometrical radius r g and the mean freeze-out time 
Tf are of course larger for PbPb than for hh collisions, but surprising is the 
similarity of the duration At/ of emission in both. The fact that At/ ~ t/ 
in hh collisions, indicates that the radiation process occurs during all the 
evolution of hadronic matter in this type of collisions. On the other hand, 
At/ < Tf for PbPb collisions suggests that there the radiation process only 
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sets in at the end of the evolution. Other important differences are the large 
transverse flow velocity (ut) and small transverse temperature gradient in 
PbPb as compared to hh collisions. 

4-6.5. The space-time distribution of ir emission: 

Figure 15a gives a reconstruction of the space-time distribution of pion 
emission points [116], expressed as a function of the cms time variable t 
and the cms longitudinal coordinate z. The momentum-integrated emission 
function along the z-axis, i.e., at (x,y) = (0,0) is given by 

It relates the parameters fitted to the NA22 data with particle produc- 
tion in space-time. Note that the coordinates (t,z), can be expressed with 
the help of the longitudinal proper-time r and space-time rapidity rj as 
(rcosh?7,rsinh77). 

One finds a structure resembling a boomerang, i.e., particle produc- 
tion takes place close to the regions of z = t and z = —t, with gradually 
decreasing probability for ever larger values of space-time rapidity. Al- 
though the mean proper-time for particle production is Tf = 1.4 fm/c, and 
the dispersion of particle production in space-time rapidity is rather small 
(A77 = 1.36), a characteristic long tail of particle emission is observed on 
both sides of the light-cone, giving more than 40 fm longitudinal extension 
in z and 20 fm/c duration of particle production in the time variable t. 

An, at first sight, similar behavior is seen in Fig. 15c for PbPb collisions 
[117]. An important quantitative difference is, however, that particle emis- 
sion starts immediately in hadron-hadron collision, but only after about 4-5 
fm/c in PbPb collisions! 

The information on («t) an d (^?) from the analysis of the transverse 
momentum distribution can be used to reconstruct the details of the trans- 
verse density profile. An exact, non-relativistic hydro solution was found 
[119] using an ideal gas equation of state. In this hydro solution, both 

AT m(u T ) 2 

{—) < — — > ( 62 ) 

are possible. The < sign corresponds to a self-similar expanding fire-ball, 
while the > sign corresponds to a self-similar expanding ring of fire (see 
Fig. 16). 



33 




Figure 15. The reconstructed emission function S(t, z) in arbitrary vertical units, 
as a function of time t and longitudinal coordinate z (left diagrams), as well as the 
reconstructed emission function S(x,y) in arbitrary vertical units, as a function of 
the transverse coordinates x and y (right pictures), for hh (upper pictures) and 
PbPb (lower pictures) collisions, respectively [116, 117, 118]. 
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Figure 16. Illustration of the development of smoke-ring solutions for large tem- 
perature gradients in exact solutions of non-relativistic hydrodynamics [119]. 

Assuming the validity of this non-relativistic solution, one can recon- 
struct the detailed shape of the transverse density profile. The result looks 
like a ring of fire in the x,y plane in hh interactions (Fig. 15b), while in 
PbPb collisions it has a Gaussian shape (Fig. 15d). 

The formation of a ring of fire in hh collisions is due to the rather small 
transverse flow and the sudden drop of the temperature in the transverse 
direction, which leads to large pressure gradients in the center and small 
pressure gradients and a density augmentation at the expanding radius of 
the fire-ring. This transverse distribution, together with the scaling longitu- 
dinal expansion, creates an elongated, tube-like source in three dimensions, 
with the density of particle production being maximal on the surface of the 
tube. 

The pion emission function S(x, y) for PbPb collisions, on the other 
hand, corresponds to the radial expansion, which is a well established phe- 
nomenon in heavy-ion collisions from low-energy to high-energy reactions. 
This transverse distribution, together with the scaling longitudinal expan- 
sion, creates a cylindrically symmetric, large and transversally homogeneous 
fireball, expanding three-dimensionally with a large mean radial component 
(ut) of hydrodynamical four- velocity. 

Because of this large difference observed for those two types of collision, 
analysis of the emission function in e + e _ collisions is of crucial importance 
for the understanding of the actual WW overlap and has been started. 

4.7. The 7r°7r° system 

In a string model, unlike Tr^vr^-pairs, pairs of prompt 7r°'s can be emitted 
in adjacent string break-ups. In momentum space, the correlation function 
is, therefore, expected to be wider for neutral pions than for charged ones. 
Neutral pions, furthermore, do not suffer from Coulomb repulsion. How- 
ever, the detection of several 7r°'s in one event requires high efficiency of 
7-detection in a wide energy range and geometrical acceptance. Further- 
more, the correlation function at small Q is strongly influenced by resonance 
decays as rj — ► 7r 7r°7r°, rj' — > 7r°7r°?7, Kg, fa — ► 7r°7r° and other final-state 
interactions [120]. 

First evidence for Bose-Einstein correlations in 7r°7r° pairs was found in 
[121]. In a first measurement of the radius in 7r~Xe interactions at 3.5 GeV 
[122], the size of the 7r° emission region was found compatible with that for 
charged pions. 

The question was taken up again by L3 [123], where both tq and A are 
found to be on the low side when compared to the tt^tt^ results obtained 
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under the same experimental conditions. The difference in A can at least 
partially be explained by the contribution of resonances. The difference in 
size parameter is r±± — vqq = 0.150 ± 0.075 ± 0.068 fm. 



4-8. Higher- Order Bose- Einstein Correlations 

4-8.1. The formalism 

It is convenient to use the normalized inclusive density and correlation 
functions already defined in Eqs.(10) and (11). The normalized inclusive 
density for two identical pions is 

i? 2 (l,2) = l + K 2 (l,2). (63) 

In the limit of a completely chaotic and static pion source, #2(1, 2) reduces 
to the square of the Fourier transform F(p 1 — p 2 , E\ — E^) of the space-time 
distribution of the source, #2(1,2) = |F(1,2)| 2 , where and (i = 1,2) 
are the three-momentum and energy of pion i, respectively. 

If the Gaussian parametrization is used for l-^XQ^)! 2 , then one has 

K 2 {Ql) = \F{Ql)\ 2 = zM-rlQl) ■ (64) 

In terms of the Qij variables and for the case of a completely chaotic 
source, the normalized inclusive three-pion density is [124, 125] 

^3(1,2,3) = iHn<&)\ 2 + \n<&)\ 2 + \n<%i)\ 2 

+ 2Re{F(Q2 2 ) F(Q 2 3)i?(Q 2 i)} j (65) 

so that the genuine three-particle correlation reads 

K 3 (l, 2, 3) = 2Re{F{Q\ 2 )F{Qlz)F(Ql l )}. (66) 

In general, the genuine three-particle correlation #3(1, 2, 3) is not ex- 
pressed completely in terms of the two-particle correlation function (64), 
but contains also new information on the phase 4>ij of the Fourier transform 
of the source, 

#3(1,2,3) 

cos© = — . , (67) 

2 V /# 2 (1,2)K 2 (2,3)# 2 (3,1) 

with (j> = 4>i2 + 4>23 + 031 being a function of Qij and cos (j) — ^ 1 as Qij — > 0. 
Geometrical asymmetry in the production mechanism (emission function) 
due to flow or resonance decays will only lead to small (few percent) re- 
duction of cos0 from unity [126]. Equation (67) is not valid for (partially) 
coherent sources and more complicated expressions are needed [126]. If 
cos cf) considerably differs from unity at Qij > 0, one can infer that partial 
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coherence is present (or, alternatively, that K3 is suppressed due to dilution 
in the case of many independent sources!). 

To the extent that phase factors may be neglected and the Gaussian 
approximation would hold, K% is related to K2 via the expression 

2 

K 3 (Q 2 3 ) = 2eM~Ql) = ^K 2 {Ql) (68) 

with 

Ql = Qi23 = (Pi + P2 + Pz? ~ 9M^ = Q\ 2 + Qj 3 + Q 2 23 . (69) 



4-8.2. Genuine three-particle correlations 

Non-zero genuine correlations up to order q = 5 were first established 
by the NA22 collaboration [127] in terms of cumulant moments. So they 
must show up here, as well. The function K^{Q^) + 1 is given in Fig. 17a. 
A non-zero K% is indeed observed in the data for < 0.2 (GeV/c) 2 [61], 
but not in FRITIOF. 
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Figure 17 a) The normalized three particle correlation function K${Q\) added to 1. 
The full line is the result of a fit by (70), the dashed line corresponds to FRITIOF 
results [61]. b) The function Kz(Q)-\-l for like-sign triplets and unlike-sign triplets. 
The predictions of JETSET without BE (dashed line) and with BE correlations 
(full line) are also shown [128]; c) Like-sign triplets after Coulomb correction, with 
a Gaussian fit (solid line) [129]. 

Both observations, the existence of genuine three-particle correlations 
and the underestimate in JETSET are supported by DELPHI [128]. In 
Fig. 17b, the three-particle correlation function K% + 1 is shown for like- 
charged triplets (upper) and unlike-charged triplets (lower), respectively, 
together with the prediction of JETSET with and without BE correlations. 
The parameters used to include the BE correlations are the same as in the 
two-particle correlation study of DELPHI [30]. The model is in reasonable 

agreement with the data for the (+ H — ) and (-1 ) configurations but 

underestimates the enhancement for the (+ + +) and ( ) correlations. 

Bose-Einstein interference in JETSET not only changes the distribution 
of like-sign correlations, but also the unlike-sign ones and leads to better 
agreement with the data. 

Statistically a better evidence for genuine three-particle BE correlations 
now comes from OPAL [129]. This is shown in Fig. 17c, together with a 
Gaussian fit over the range 0.25 < Q3 < 2.0 GeV, giving r 3 = 0.580 ± 
0.004 ± 0.025 fm and A 3 = 0.504 ± 0.010 ± 0.041. Within two standard 
deviations, the value for r 3 agrees with the relation r 3 = r 2 /V2 (see (68)) 
when compared to r 2 obtained in [77]. 

The question is, whether the observed genuine three-particle correlation 
can be fully expressed in terms of the simple product of two-particle corre- 
lation functions according to (68), or whether information can be extracted 
on the relative phases of (66). If relation (68) holds, the function l+K%(Q\) 
can be described by the parameters r 2 = 0.85 ±0.01 fm and A2 = 0.38 ±0.02 
deduced from the fit of the normalized two-particle density R2{Q 2 ): 

K 3 (Ql) + l = 7 [l + 2Xl /2 exp(-^r 2 2 Q 2 3 )](l + 5Ql) . (70) 

Within the errors of NA22, the resulting parameters r 2 and A2 do not con- 
tradict those of the two-particle correlations and, therefore, are consistent 
with Eq.(67) and, therefore, with incoherent production of pions. DELPHI 
and OPAL unfortunately did not make use of this possibility, but L3 did 
[131]: 

Fig. 18 gives cos (ft (Eq.(67)) as a function of Q3 for the case that the 
cumulants K 2 and K3 are parametrized in terms of a first-order Edgeworth 
expansion of a Gaussian. The L3 result is consistent with cos <p = 1 for all 
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Qs and therefore with full incoherence. 
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Figure 18. cos0 as a function of Q3 assuming R2 is described by the first-order 
Edgeworth expansion of the Gaussian [131]. 

Three-pion correlations have also been studied in heavy-ion collisions 
[132, 133] and (cos^) = 0.20 ± 0.02 ± 0.19, i.e. no genuine three-particle 
correlations are found outside the (large) errors for SPb [133]. The authors 
interprete this result as evidence for partial coherence [133]. 

What is particularly remarkable, however, is that the same experiment 
(NA44) with the same methodology finds an average {cos (ft) = 0.85 ± 0.02 ± 
0.21 for PbPb collisions [133] and that this is supported by a value of 
(cos0) = 0.606 ± 0.005 ± 0.179 earlier reported by WA98 [132]. 

So, if we trust NA44 (and I have no reason not to) and try to stick with 
conventional pion interferometry, we end with a beautiful dilemma: 

i) e + e~ collisions are consistent with fully incoherent pion production 
(cos 4> ~ 1)! 

ii) SPb collisions are consistent with coherent production (cos0 ~ 0)! 

iii) PbPb is somewhere in between! 

It could not be more opposite to any reasonable expectation from con- 
ventional interferometry [134]. The hint for an alternative interpretation 
comes from a comparison of Eqs. (67) and (46). What conventional inter- 
ferometry calls the cosine of a phase has in fact nothing to do with a phase. 

3 /2 

It is simply the ratio of K% and twice K 2 ■ It may be a challenge for the 
string model to explain why this is unity for an e + e~ string. If that can 
be explained, the rest looks easy and very much in line with the behavior 
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of the strength parameter A discussed at the end of Sect. 4.5: The ratio 
cos <fi ~ K^/2K 2 decreases with the number of independent sources N like 
N 2 /2N 3 / 2 oc N 1 / 2 . As A does, it decreases with increasing atomic mass 
number A up to SPb collisions. The saturation or increase of A at and 
above this A has been explained by percolation [109] of strings in Sect. 4.5. 
Exactly the same explanation can be used to understand an increase of the 
ratio (not the phase!) cos 4> between SPb and PbPb collisions. 

5. Conclusions 

In view of possible inter- W Bose-Einstein Correlations distorting fully 
hadronic WW final states in e + e~ collisions, the state of the art has been 
summarized on Bose-Einstein correlations in Z fragmentation. Where not 
(yet) available from the Z, information has been borrowed from other types 
of reactions. We consider this experimental information a major challenge 
to existing and future models. 

1. Bose-Einstein correlations definitely exist in Z fragmentation. There 
is no reason that they should not exist within a single W (intra-W BEC). To 
understand possible presence or absence of BEC between pions originating 
from different Ws (inter-W BEC), profound knowledge of BEC in the high 
statistics Z fragmentation data is obligatory. 

2. The correlation is far from spherically symmetric. The correlation 
domain (defined by the lengths of homogeneity in three space directions) is 
elongated along the event axis. Because of strong space-momentum corre- 
lations, this elongation is small, however, as compared to the length of the 
total string. 

3. The correlation is far from Gaussian. Good results have been obtained 
with an Edgeworth expansion, but even power-law behavior is not excluded. 

4. A scaling first observed in heavy-ion collisions is now also 
observed in Z fragmentation and may suggest a "transverse flow" even there. 

5. Dilution of correlations in high- multiplicity hh and A A collisions 
suggests the lack of cross talk between neighboring strings at low density of 
strings, but percolation may set in at the highest densities. 

6. The full emission function in space-time can be extracted from a com- 
bination of inclusive single-particle distributions and BE correlation func- 
tions. So far, this has only be done in hh and AA collisions, in the framework 
of a model for three-dimensional hydrodynamic expansion. While a Gaus- 
sian shaped fireball is observed for AA collisions, a firetube is observed for 
hh collisions. A study of e + e~ collisions is under way. 

7. Consistently with the expectation from the string model, the radius 
for 7r°7T correlations is found to be smaller than that for 7r ± 7r ± correlations. 
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8. Genuine three-particle correlations exist in Z fragmentation and, 
according to conventional interpretation, would allow to extract a "phase" 
unmeasurable in two-particle correlations. The resulting zero phase in e + e~ 
is consistent with what would be expected for fully incoherent emission. 
Comparison to the results obtained for heavy- ion collisions, however, raises 
doubts on the conventional interpretation. 
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